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Pressure anisotropy and small spatial scales induced by velocity shear 
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Non-Maxwellian metaequilibria can exist in low-collisionality plasmas as evidenced by satellite 
and laboratory measurements. By inclnding the full pressure tensor dynamics in a fluid plasma 
model, we show that a sheared velocity field can provide an effective mechanism that makes an 
initial isotropic state anisotropic and agyrotropic. We discuss how the propagation of magneto¬ 
elastic waves can affect the pressure tensor anisotropization and its spatial filamentation which are 
due to the action of both the magnetic field and flow strain tensor. We support this analysis by a 
numerical integration of the nonlinear equations describing the pressure tensor evolntion. 


The aim of this Letter is to show that a sheared velocity 
field in a weakly collisional, magnetized plasma drives a 
macroscopic pressure anisotropization in the plane of the 
velocity strain tensor. This represents a general mech¬ 
anism when collisional relaxation is absent or slow that 
causes part of the kinetic energy of the plasma flow to 
be locally transformed into anisotropic “internal energy”. 
This energy conversion implies that shear flows do not af¬ 
fect the plasma dynamics only through the fluid destabi¬ 
lization of Kelvin-Helmholtz (KH) modes [I| or by break¬ 
ing the correlation length of unstable modes [1,3 respon¬ 
sible e.g., for anomalous energy transport in magnetically 
confined plasmas, but can lead to the onset of additional 
phase space instabilities driven by the induced pressure 
anisotropy. 

In magnetized plasmas the fast particle gyromotion 
in a sufficiently strong field makes the pressure tensor 
isotropic in the plane perpendicular to the magnetic field 
direction but allows for different parallel and perpendic¬ 
ular pressures (gyrotropic pressure as is the case for the 
double-adiabatic or CGLQ closure). On the contrary, 
the fluid strain = dui/dxj in the sheared fluid veloc¬ 
ity Ui{x) has a twofold effect: first, through its rotational 
component it combines or competes with the gyrotropic 
effect due to the magnetic field, second it induces pres¬ 
sure anisotropy (agyrotropic pressure) in the plane per¬ 
pendicular to the magnetic field (taken to coincide with 
the velocity shear plane) through its incompressible rate 
of shear (its symmetrical traceless component). 

Here we discuss the role of the flow strain in the dy¬ 
namic equations of the full pressure tensor as obtained 
from the second moment of Vlasov Equation (VE), thus 
going beyond both the CGL closure and the Finite Lar- 
mor Radius (FLR) corrections 3 approach. We focus in 
particular on the dynamics of the full pressure tensor 
within a 1-fluid description of a dissipationless magne¬ 
tized plasma and show how the propagation of “magneto- 
elastic" waves can affect the pressure anisotropization 
and small spatial scale formation due to the interplay 
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between the gyrotropic and the non-gyrotropic dynamics 
induced by the magnetic field and by the strain tensor. 

Non-Maxwellian states, sometimes exhibiting pressure 
agyrotropy^-Q, are observed both experimentally [6l-[l3| 
and in Vlasov simulations [Tsl IT^ ■ the role of a shear 
flow in affecting the kinetic properties of a collisionless or 
weakly collisional plasma is well established experimen¬ 
tally. Even if mechanisms based on the GGL paradigm 
are now almost acquired when explaining the main 
features of solar wind anisotropies jl^ . the correlation 
between the presence of a velocity shear and the extent 
of anisotropy in particle distributions is evoked both 
for the core protons in the fast solar windl^ and 
in “space simulation” laboratory experiments [l4| . A 
sheared convection velocity in the Earth ionosphere is 
thought to play a role in the heating of ions and in the 
consequent plasma bulk upflow in the auroral region (T^. 
Moreover, the presence of a velocity shear is known to 
play an important role in the enhancement of a variety 
of pressure anisotropy-related plasma instabilities. The 
presence of a velocity shear in the near-Earth plasma 
sheet profile prior to a substorm expansion lowers 
the instability threshold of ion-Weibel modes in the 
geomagnetic tailfl^. The role of a vorticity-related 
velocity shear in generating a gyrotropic temperature 
anisotropy in ion temperature gradient (ITG) driven 
turbulence, was discussed in Ref. 3 in connection 
with the relaxation of transport barriers in tokamak 
plasmas. Finally, an anisotropic pressure tensor was 
shown to form during the nonlinear stage of the current- 
filamentation instability (GFI) arising in the presence of 
two opposite cold electron beams. This anisotropy was 
shown to increase the threshold and growth rate of the 
reconnection instability developing on the shoulder of 
the GFI-generated magnetic structuresfl^. These also 
develop in the presence of radially inhomogeneous beams 
such as in high intensity laser-plasma interactions [20l| 
and are measured in laboratory experiments . 

In Ref. [ 2 ^ sufficient conditions for the instability 
of electromagnetic waves in an unmagnetized Vlasov 
plasma with a sheared velocity distribution were studied, 
showing that an anisotropic pressure tensor leads to 
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the instability of transverse perturbations. This was 
applied in to study the steady state achievable when 
the shear-induced temperature anisotropy balances 
the electron diffusion in the velocity space due to the 
growth of static magnetic perturbations generated by 
the pressure anisotropy itself. Anisotropic turbulence 
induced by a KH unstable velocity shear and by a 
Von Karman flow[2^ was also pointed out. 

We start from the two-fluid equations of a colli¬ 
sionless magnetized plasma, obtained by evaluating the 
moments of VE. We neglect the electron dynamics (i.e. 
rue/mi —?► 0) and temperature, whereas the components 
of the full pressure tensor contribute to ion dynamics. 
The 2”'^ anisotropic ion moment, 11^^, is closed by 

= — £u(n^j) -f A^u(nij) (1) 

jVu|=T“^ 

where = q\B\/{me) and we have introduced the linear 
operators = d]^{u}^Y\.j^j) -t-and 

SjijyiTl{iBjyi)j and the cor¬ 
responding characteristic hydrodynamic {t^^ = |Vu|~^) 
and magnetic (r^ = time scales. The simplifying 
assumption, frequently used in the literature , of 

neglecting in Eq.(IT]) the divergence of the ion heat flux 
tensor, Qijk = {mn{vi — Ui){vj — Uj){vk — Uk)), is consis¬ 
tent with the geometrical configuration considered later 
in this letter at least until very short spatial scales in 
the plane perpendicular to the magnetic field are formed 
during the nonlinear evolution. Typical closures of II 
have been performed by means of a power expansion 
each in terms of some small parameter: including cor¬ 
rections in Eq. o due to a small collision time with re¬ 
spect to both Tjj and , leads to Braginskii’s gyrovis- 
cous FLR modelfl^, while a small leads to FLR 

gyrotropic corrections to CGL equations (sl WA. . More 
recent FLR-Landau-fluid models [13 also retain Landau- 
fluid effects [3l|. Here we consider a dissipationless regime 
and do not assume the ratio to be small. 

Defining the matrices B^- = SijmBm, and Wy = 
{diUj — djUi)j2 that describe the rotation induced by 
the magnetic field and by the shear flow respectively, the 
strain traceless matrix = {djUi -I- diUj)f2 — CSij, the 
compression C = {dkUk)/?> and the derivative d/dt = 
dt + Ukdk, Eq. m can be conveniently written as 

^n = [B-hW,n]-{D,n}-5C'n ( 2 ) 

where [, ] denotes commutator and { , } anticommutator. 
The first r.h.s. term shows that the magnetic field B and 
the flow vorticity u: (wi = SijkWjk) combine to make II 
to rotate around the axis of B -|- W. The perpendicu¬ 
lar components rotate at twice the cyclotron frequency 
in absence of vorticity, or at twice the fluid rotation fre¬ 
quency in the unmagnetized case. If the axes of B and 
W are aligned the two frequencies add up if B • a; > 0 


and subtract if B • u) < 0. The last r.h.s. term of Eqs. © 
acts isotropically on II while the second term can induce 
pressure anisotropization whenever D is not zero. 

First, we consider a model system with an incompress¬ 
ible shear flow Uy{x) constant in time (energy is thus 
constantly injected from outside) in the presence of a 
uniform and constant magnetic field along the z-axis. In 
this model the velocity strain and the vorticity have the 
same magnitude, B is uniform in space, the axes of B 
and W are aligned along z while D has no z components. 
Eq.(I2]) leads to a linear system in which the components 
of n perpendicular to B are decoupled from those hav¬ 
ing a parallel component and their determinant is con¬ 
stant. Three eigenvalues are obtained for the perpendic¬ 
ular components: 70 = 0 , which corresponds to a station¬ 
ary mode with and Il].y = 0 , and 

7 ± = ±2i^Jflcfl'{x) with Hyy/n 2 x = —Q!{x)/Q.c and 
BZt= ±* ^J^l'{x)/^lc■ Here H'(x) = + dxUy{x). 

Provided ^l’{x) > 0, the 70 mode can describe an 
equilibrium solution of Eq.Q (in agreement with the 
self-consistent equilibria discussed in [13), Byy/B^x = 
H'(x)/Hc and = 0. The 7 ± modes represent either 
oscillations or growing and damped modes depending 
on the sign of Q,'{x). For H'(a;) > 0 the perpendic¬ 
ular pressure tensor components of an initial isotropic 
state with na;a;(a;,0) = nyy(a;,0) = P±{x) oscillate 
in time around a mean value given by (Jiyy{x^t)) = 
{Vt'{x)/Vtc){Bxx{x,t)) = {Q.'{x) + Q.c)Px{x)/{2Vtc) and 
lJVxy{x:t)) = 0, the amplitude of the oscillations of 
Iiyy{x,t) being dxyPy{x)P^_{x,0)/{2Vl{x)). In FigUlthe 
profile of Ilyy is shown at different times, for an initial 
pressure tensor Hy = Sij, = 1 and Uy = Vocos{kx) 
with Vo = 1-5 and k = 1. An important feature caused by 
the spatial inhomogeneity of the shear flow is the strongly 
inhomogeneous growth of the components of the pres¬ 
sure tensor, as regions where the evolution is oscillatory 
alternate, depending from the local sign of il'il, with re¬ 
gions of exponential growth occurring over a time scale 
T„ = {kVo)~^. This gives rise to a spatially filamented 
pressure tensor. 

Second, we consider the self-consistent (SC) case in 
which the flow and the e.m. fields evolve in time 
according to Eq. o and to the ideal MHD equations 
with an anisotropic ion pressure (note the Hall term in 
Ohm’s law to vanish identically for the examples dis¬ 
cussed in this Letter). This system conserves the to¬ 
tal energy (pu^/2-|--|-tr{n}/2), 

and depends on three dimensionless parameters th/tb 
= {cAlcH){LH/di), {c^jc^f and {cBjc^f with the 
scalelength of the configuration, ca the Alfven velocity, 
= Puhn 8 - measure of the flow velocity, di = c^/ftc 
the ion skin depth, and c± = P±/p = Csf'/2 with Ce 
the “sound” velocity evaluated with respect to the initial 
ion pressure, assumed isotropic in the plane perpendicu¬ 
lar to B [33|. Two parameters only, {t^jIt^) j{c^jcjj)‘^ 
and {tjj jxj^)f{c±f , rule the linear dynamics. The 
CGL-FLR limit of Ref.Q is recovered in the low fre¬ 
quency limit dt He for {c^/Cj^){di/L <C 1 and 
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dijLjj ^ {c^jCj^dijL^ 1, where no specific order¬ 
ing for cx/c _4 is assumed. 

In the SC case the anisotropization of the pressure 
tensor caused by the presence of an initially imposed 
shear flow is limited by the the reaction of the pressure 
tensor on the plasma flow which reduces its shear and 
by the excitation of nonlinear ‘‘‘‘magneto-elastic'' pertur¬ 
bations that tend to propagate the shear of the veloc¬ 
ity flow outwards. The main features of these pressure 
tensor and velocity perturbations can be understood by 
referring to the linear waves described by the SC sys¬ 
tem that propagate in a uniform, homogeneous equilib¬ 
rium with B = BqBz, density poi double adiabatic pres¬ 
sure Bij = P±6ij (P|| — P±)BiBj/B‘^, and wave-vector 
k = kex- The resulting dispersion relation 

+ 3c^)]—= 0, 

( 3 ) 

has a higher frequency branch HFB, uif^ ~ AB,^-\-2k^c\ 
and ^ +3c5^) for k'^ —>• 0, and —>■ oo re¬ 

spectively, and a lower frequency branch LHB, uif ^ 
k‘^{c\ 2c']_) and ujf ~ k'^cj_ for —>• 0, and k'^ —>• oo. 

The —>■ 0 limit of the LFB corresponds to the CGL 
form of a perpendicular magnetosonic wave. In the limit 
of vanishing magnetic field the two branches are not dis¬ 
persive and reduce to a longitudinal and to a transverse 
sound mode propagating at phase velocities v^cj, and 
c_L respectively. Note that in this latter limit the SC 
system admits the propagation of finite-amplitude trans¬ 
verse waves where Uy and 11^^ satisfy the wave equa¬ 
tion df = c]_d^ and the energy continuity equation 
dt{poUy Byy) = —dx(uyByy). Whon excited, these 
waves may carry away and disperse an initially imposed 
shear flow Uy{x) with spatial scale on time scales 
of the order of ~ thus reducing the 

anisotropization and the spatial inhomogeneity of the 
pressure tensor forced by the flow. This propagation 
might affect the onset of the KH instabilities in a colli¬ 
sionless weakly magnetized plasma. In the presence of a 
strong magnetic held these waves, in particular the HFB, 
couple the x and the y components of the plasma veloc¬ 
ity. In addition they become dispersive with the group 
velocity of the HFB going to zero for fc —>■ 0. 

The nonlinear SC case has been integrated numeri¬ 
cally starting from an isotropic initial condition with 
homogeneous density, B = B^e^ and u = M°(x)ey, 
varying the value of the ratios of the three dimen¬ 
sionless parameters. In Figsl^IS] we consider the case 
with Uy{x) = Vb tanh(a;/di)/cosh^(a;/di) and = 

^aI^h = = 1- The results obtained can be qual¬ 

itatively accounted for by referring to the linear modes 
described above where the initial shear velocity Uy{x) is 
interpreted as an initial perturbation. Note that, though 
its characteristic scale-length is chosen of the order di , the 
initial Fourier spectrum peaks around kdi < 1 (Figl2|). 
The initial perturbation can be written as a superpo¬ 
sition of the LFB and of the HFB. To leading order 
in kdi 1 the polarization vectors components in the 


(ux,Uy) basis are {1,—i o{kdi)) and (I,—i) for w; and 
for LOh respectively. This implies that the chosen ini¬ 
tial perturbation corresponds to a superposition of the 
two branches with equal and opposite amplitudes and 
that the time evolution of Uy{x) is mainly determined 
by that of the HFB. This is consistent with the re¬ 
sults of the numerical integration and explains why the 
pressure agyrotropy (Figs 13]), that in our geometry is 
mainly related to the spatial inhomogeneity of Uy{x), 
tends to remain at the original position and not to be 
carried away at the Alfvenic group velocity of the LFB, 
at least until small spatial scales are formed which are 
instead transported away efficiently by the HFB. For ex¬ 
ample, at a; = 0 Figs.® show local peak agyrotropies 
\Byy — Bxx\/\B^y-\-B°xx\ of 0.12, 0.22 and 0.14, and “mean 
anisotropies” 2Bzz/(Byy + Bxx) of 1.02, 2.35 and 0.98 for 
the = 0.1,1,10 case respectively. On the contrary, 

both branches contribute to the evolution of Ux (x) where 
the initial cancellation is removed as time evolves with 
the LFB component propagating outwards and the HFC 
essentially mirroring, with an inverted sign, the behavior 
of Uy . An increase of the magnetic field Bq has a double 
role: on the one hand it tends to enforce perpendicular 
gyrotropy while on the other, when the ratio c^/c_l is 
increased, the group velocity of the HFB decreases and 
the initial perturbation of Uy remains longer confined at 
its initial position. 

The fluctuations of Ux being compressible induce fluctua¬ 
tions of Bz and Il^z consistent with the magnetoacoustic 
polarization bH^^/H^^^ = SB^/B^ (non shown here). 

The interplay between the filamentation shown in Fig]!] 
and the propagation of the disturbances of the pressure 
tensor result in the formation of fine-scale spatial struc¬ 
tures. The latter dominate for higher values of Bq be¬ 
cause of the nonlinear steepening of the front of the prop¬ 
agating “magneto-elastic waves”, as evidenced by the Ux 
and Ilxx profiles in Figs|3||3| 

The following points appear as particularly relevant in 
the analysis of the SC case: 

1) analogously to the CGL paradigm, the ratio between 

the average pressure in the perpendicular plane and the 
parallel component can vary considerably (e.g. = 

1 case in FiglSj), 

2) agyrotropic anisotropization also occurs, though less 
pronounced as the magnetic held increases, since the flow 
energy may be preferentially transferred to one tensor 
component only in the plane perpendicular to the 
magnetic held direction. 

3) The pressure anisotropization is spatially asymmetric, 
depending on the sign of cu • 15, and tends to be localized 
in space while moving together with the initial velocity 
inhomogeneity which is carried away by the nonlinear 
magneto-elastic type wave-packets, 

4) a relaxation toward new non-MHD equilibria like 
those discussed in may be asymptotically achieved, 
even if, at the increase of Bq, the spatial inhomogeneity 
u = Uy{x)ey transferred to the relevant Hy components 
takes longer to leave the initially sheared region when it 
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FIG. 1. Evolution of \lyy{x,t) for B = BqSz and constant 
u — (0, Vb cos(a;/cii), 0), Q.cT^ = 1 and Vo = 1.5c^. Both the 
exponential growth > 0) and the spatial filamentation 

of the oscillating solutions {^'{x) < 0) are visible. 


is peaked at large wavelengths {kdi < 1). 

The plasma dynamics described above explains why 
isotropic MHD equilibria cease to be equilibria in pres¬ 
ence of a stationary sheared flow and why an initial 

anisotropic distribution function is needed to initialize ki¬ 
netic simulations in presence of a velocity shear [s^. This 
can affect the onset and development of the KH instabil¬ 
ity and may have implications for the evolution of trans¬ 
port barriers in tokamak turbulence 0. Another direct 


implication is for turbulence itself: since small-scale spa¬ 
tial inhomogeneities are naturally developed during the 
direct cascade, we may expect that isotropic turbulent 
states are not likely to exist whenever a full pressure ten¬ 
sor evolution is accounted for. In particular, since non- 
negligible discrepancies with respect to the CGL closure 
become important when ~ flcj for ^ (Alfvenic 
turbulence) pressure anisotropies in the plane perpendic¬ 
ular to the magnetic field can be expected when velocity 
inhomogeneities are generated at a scale ^ di, appar¬ 
ently in agreement with the temperature anisotropiza- 
tion observed in Refs.[l5l| in concomitance with the de¬ 
velopment of current and vorticity layers of thickness 
~ di- Finally we recall that the occurrence of an agy- 
rotropic pressure tensor is well documented in solar wind 
measurements!^ 01, possibly correlated to plasma flows, 
see e.g. Ref.@. 
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10 respectively (here = di). 


















